Induced superfluidity of imbalanced Fermi gases near unitarity 
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The induced intraspecies interactions among the majority species, mediated by the minority 
species, is computed for a population-imbalanced two-component Fermi gas. Altliough the Feshbach- 
resonance mediated interspecies interaction is dominant for equal populations, leading to singlet 
s-wave pairing, we find that in the strongly imbalanced regime the induced intraspecies interaction 
leads to p-wave pairing and superfluidity of the majority species. Thus, we predict that the observed 
spin-polaron Fermi liquid state in this regime is unstable to p-wave superfluidity, in accordance with 
the results of Kohn and Luttinger, below a temperature that, near unitarity, we find to be within 
current experimental capabilities. Possible experimental signatures of the p-wave state using radio- 
frequency spectroscopy as well as density-density correlations after free expansion are presented. 

PACS numbers: 05.30.Fk, 03.75.Ss, 67.85.-d, 32.30.Bv 



I. INTRODUCTION 

The extraordinarity variety of tunable "knobs" in cold 
atomic gas experiments has yielded a wide assortment 
of correlated phases of matter ranging from superfluid 
and Mott insulating phases of bosonie atoms, to super- 
fluid and nonsuperfluid phases of fermionie atoms [1-4]. 
In the fermionie case the experimentally available knobs 
include the interactions among two species of fermion, 
their relative densities [5-8], and also the effective spa- 
tial dimension, which can be controlled via an appropri- 
ate confining trapping potential. 

Our present focus is on the case of a three-dimensional 
two-species Fermi gas (labeled by the spin index a =t, i) 
interacting via a magnetic field-tuned Feshbach reso- 
nance. In the balanced case of equal densities of the two 
species, and as a function of the Feshbach resonance de- 
tuning that essentially controls — with a^ the s-wave 
scattering length, the atomic gas undergoes the well- 
known crossover from a Bose-Einstein condensate (BEC) 
of tightly bound molecular pairs at -r^ — > 1 , through 
the strongly correlated unitary regime at -^^ ~ 0, to a 
weakly coupled Bardeen-Cooper-Schrieffer (BCS) super- 
conductor of Cooper pairs at jT^ ^ — 1- Here, fcp oc n^'^ 
is the Fermi wavevector (with n the total density), the 
inverse of which characterizes the typical interparticle 
spacing. Importantly, the low-temperature state of this 
balanced Fermi gas is believed to undergo no symmetry- 
changing phase transitions as the BEC-BCS crossover is 
traversed. 

The behavior of strongly interacting 3D Fermi gases in 
the imbalanced case is considerably richer, with numerous 
phases having been predicted [9-16] to occur (as recently 
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reviewed in Refs. 17, 18), including the long sought-after 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state [19, 20], 
in which the superfluid spontaneously develops a spatial 
modulation in the local pairing amplitude to accomodate 
the excess fermions of the majority species [21-23]. Un- 
fortunately, the FFLO phase has not yet been observed 
in 3D imbalanced Fermi gases (although evidence for a 
ID analogue of this state has been found [24]). Instead, 
under an imposed density imbalance, characterized by 
the population imbalance P = ^^^^^ (note we always 

assume the densities ria- of species a satisfy n-|- > nj^), in- 
teracting Fermi gases enter a regime of phase separation 
at moderate P followed by an imbalanced nonsuperfluid 
phase at large P, as seen in the phase diagram Fig. 1. 

In the strongly imbalanced limit (P — !■ 1^) our system 
amounts to considering the phase resulting from adding a 
few spins-4. to a Fermi sea of the spins-f. In this strongly 
imbalanced regime, within the simplest picture there are 
two possible fates of such an added spin down: It can 
cither form a molecular bound state with one spin-^, or 
it can remain unpaired. In the former case, occurring at 
fls > but kpas ^ 1 (i.e. the BEC limit), as more spins- J, 
are added more molecules will develop and, presumably, 
condense. Our interest is in the latter regime, in which 
(within this simple picture) no interspecies pairing occurs 
and condensation is suppressed. 

From one point of view, the inability to establish pair- 
ing in this regime is simply traced to the fact that, with 
an imposed density imbalance, not all the majority spins- 
'[ have a minority spin- J, fermion to pair with. More 
precisely, the presence of a density imbalance implies a 
concomitant Fermi surface (and Fermi energy) mismatch, 
so that the formation of low-energy Fermi surface pairing 
correlations is interrupted, leaving an imbalanced inter- 
acting Fermi liquid state [25-27] that, apparently, lacks 
any superfluid order even at low temperatures. 

How do the strong attractive interspecies interactions 
manifest themselves, given the inability of this system to 
form s-wave Cooper pairs and condense? Recent experi- 



ments show evidence of the formation of Fermi polarons, 
in which a cloud of the spins-t form around each spin- 
J,, leading to an observable shift in the spin- J, chemical 
potential [28, 29] in agreement with theoretical predic- 
tions [30-32] . The central question studied here concerns 
whether, at T — > 0, this polaron phase of matter persists 
or whether another broken-symmetry phase emerges. 

One motivation for this possibility is the well-known 
results of Kohn and Luttinger [33, 34], who showed that 
interacting Fermi liquid phases are generally unstable to 
pairing in some angular momentum channel at low T. 
Given that the polaron state of imbalanced Fermi gases 
is, at its heart, essentially a Fermi liquid for both species 
(exhibiting, for example, a sharp Migdal discontinuity 
in the momentum occupation at the Fermi surface), we 
generally expect the Kohn-Luttinger mechanism to hold, 
yielding odd angular momentum intraspecies Cooper 
pairing at both the spin-t and spin-^ Fermi surfaces for 
T ^ [35, 36]. The odd-angular momentum pairing 
requirement follows from the Pauli exclusion principle; 
while the simplest such state has p-wave symmetry (as- 
sumed here), more generally any odd £ is possible. 

In fact, p-wavc pairing is also possible in the deep BEC 
limit of imbalanced gases via a somewhat different mech- 
anism [35, 37]. As noted above, in the deep BEC limit 
tightly bound molecular pairs form even at very large 
imbalance (only vanishing when n^ — > 0). If we imagine 
decreasing P from unity in this limit, within a mean- 
field picture two possibilities emerge [12, 14, 15]: Firstly, 
the system can form a homogeneous polarized magnetic 
superfluid phase (SFm in Fig. 1), in which such molecu- 
lar pairs coexist with a Fermi sea of the excess spins-f. 
This possibility is found in the very deep BEC limit (i.e., 
■jr^ — > 2.37 within mean-field theory, with the tricritical 
point of Fig. 1 becoming a quantum tricritical point [14] 
at j^ = 2.37). Secondly, closer to unitarity, the system 
can form a phase separated mixture of SFm and imbal- 
anced normal phase. Our point here is that, in the SFm 
phase, p-wave pairing can be induced among the excess 
spins-t by the molecular bosons, a mechanism studied 
in Refs. [35, 37]. Here, we work away from the regimes 
of phase separation and SFm, focusing on induced pair- 
ing of the spins-f mediated purely by the presence of the 
spins- J, (i.e., without any molecular pairing); however, 
the relationship between these two regimes is an inter- 
esting problem for future research. 

In this paper we present the details of our calculation 
of the transition temperature below which p-wave pairing 
is expected to occur within this mechanism, expanding 
upon our recent Rapid Communication [36]. This calcu- 
lation can be summarized by the phase diagram Fig. 1 for 
imbalanced Fermi gases at unitarity, showing regions of 
polaron Fermi liquid, imbalanced or "magnetic" super- 
fluid SFm, and phase separation. Here, the first-order 
phase boundary enclosing the regime of phase separation 
and the second-order phase boundary separating the SFm 
and polaron Fermi liquid phases are only sketched (i.e. 
are not the result of a calculation). However, we drew 
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FIG. 1: (Color online) Proposed phase diagram of popula- 
tion imbalanced Fermi gases, at unitarity ^ ^^ = 0, as a 
function of temperature T (normalized to the Fermi energy 
ep multiplied by the Boltzmann constant fce) and popula- 
tion imbalance P, showing regions of imbalanced superfluid 
SFm, phase separation (at temperatures below the tricritical 
point, a solid dot (green online) [14, 42, 43]) polaron Fermi 
liquid, and a p-wave paired state. The principal question is 
whether the polaron Fermi liquid phase, in which both spins-f 
and spins-^ possess Fermi surfaces (despite the strong interac- 
tions), persist to T ^ 0. The black phase boundaries are only 
sketched but consistent with the experimentally determined 
phase diagram of Ref. 38. Below the blue curve, derived here, 
we find an instability towards p-wave pairing of the spin-f 
Fermi surface. 



these phase boundaries to be consistent with the exper- 
imental results of Ref. 38 to illustrate the fact that our 
maximum predicted transition temperature for p-wave 
pairing (solid curve, blue online) is not much smaller 
than the temperature scales characterizing these phase 
boundaries, suggesting that it may be possible to observe 
p-wave pairing of strongly imbalanced Feshbach-resonant 
Fermi gases. 

We remark that other recent work has considered the 
effect of intraspecies interactions in imbalanced Fermi 
gases [39]; however, this work assumed an intrinsic in- 
traspecies interaction (i.e., in the Hamiltonian), while our 
work assumes a vanishing bare intraspecies interactions. 
The present problem of induced intraspecies interactions 
in imbalanced gases has previously been studied by Bul- 
gac and collaborators [35], and by Nishida, the latter in 
the two-dimensional limit [40]. Additional recent work 
has studied the problem of polaron-polaron interactions 
in imbalanced Fermi gases [41], although this work did 
not address the possibility of p-wave pairing. 

This paper is organized as follows. In Sec. II we start 
from the standard one-channel model Hamiltonian for 
imbalanced Fermi gases and develop a general Green's- 
function formalism for studying intraspecies pairing in 
this setting. In this section we show that intraspecies 
pairing among the spins-f is reflected by an off-diagonal 
component of an associated Nambu self energy that. 



in turn, depends on a matrix vertex function that we 
proceed to approximate in the subsequent sections. In 
Sec. Ill, we approximate this vertex function to quadratic 
order in perturbation theory; our results in this section 
are consistent with the original calculations of Kohn and 
Luttinger [33, 34] (who studied the general problem of in- 
duced interactions in fcrmions) and, more recently, Bul- 
gac et al. [35] (who studied induced interactions in the 
present setting of imbalanced Fermi gases). Importantly, 
the perturbative formula for the transition temperature 
is invalid in the unitary regime where cold-atom experi- 
ments focus; in the regime where it applies, this transi- 
tion temperature is orders of magnitude too small to be 
observable. In Sec. IV, we attempt to go beyond the per- 
turbative limit to estimate the transition temperature for 
p-wavc pairing among the majority species in the unitary 
regime. The classes of diagrams we keep for the vertex 
function include both ladder and crossed ladder diagrams 
(which, as we show, recover the conventional Fermi liquid 
theory of imbalanced gases) , but also diagrams contain- 
ing both ladder and crossed-ladder subdiagrams (the lat- 
ter required to have a nonzero transition temperature). 
In Sec. V, we use our results from Sec. IV to estimate 
the magnitude of the pairing gap at the spin-f Fermi 
surface at T ^> 0. In Sec. VI, we turn to the question of 
how the p-wavc phase of imbalanced Fermi gases could 
be observed in radio-frequency spectroscopy and density 
correlation experiments before concluding in Sec. VII. 



II. MODEL HAMILTONIAN AND FORMALISM 

Our aim is to derive the effective interactions among 
one species of fermion, mediated by the other species 
of fermion, in a strongly imbalanced Fermi gas. In the 
present section we begin by developing a Green's function 
formalism to address this problem. Our starting point is 
the following one-channel model Hamiltonian for a gas of 
two species of fermions (labeled by a =t, |) interacting 
via an s-wave Feshbach resonance [2] (note we take h, = 



1): 



H = J2jd'r^Ur) 



2m 



*.(r) 



+ A / d^rn^{r)hi{r), 



(1) 



where the density fi„{r) = ^jj(r)^^(r) and A is the 
strength of a short-ranged pseudo-potential (approxi- 
mated by a delta- function in real space). Here, /io- is 
the chemical potential of species a; the density imbal- 
ance can be considered to arise from an imbalance in the 
chemical potentials fi^ 7^ Mi- 

The Hamiltonian Eq. (1) must be defined along with 
a cutoff reflecting the short-distance properties of the 
real physical interaction, with a corresponding scale d. 
Equivalcntly, this problem has a large momentum cutoff 
A « 2Tr/d that regularizes any divergent ultraviolet (UV) 



behavior coming from the singular nature of the delta- 
like pseudo-potential. In practice, as is well known [2], 
this can be handled by exchanging the bare coupling A 
for the vacuum scattering length Og that are related via 



47ra, 



^ 1 

T—- 



(2) 



In the weak-coupling BCS limit A — >■ 0—, this equation is 
solved by A = iis^^ j e, we can neglect the final term on 
the right side of Eq. (2). In the unitary regime where ag 
becomes large, we use a different procedure: If we assume 
physical observables are independent of A, then to study 
systems at fixed values of a^ (which is experimentally 
controllable) it is valid to replace A by a^ [using Eq. (2)] 
in approximate theoretical expressions and then take the 
limit A — >■ cxD . This strategy, which is equivalent to taking 
A -^ 0— and h. ^- 00 while holding Cg fixed via Eq. (2), 
will naturally lead to our inclusion of certain classes of 
Feynman diagrams. 

As we have discussed, the Kohn-Luttinger result im- 
plies that both the spin-^ and spin-j, Fermi surfaces of 
an imbalanced Fermi gas arc unstable at T — > 0. This 
low-T regime, possibly consisting of two interpenetrating 
p-wave superfluids, is beyond the scope of this manuscript 
and we leave it for future work. Instead, we assume that, 
owing to the large population imbalance, the transition 
temperature for any p-wave pairing of the minority spins- 
\r is much smaller than the corresponding temperature for 
the majority spins-f, as previously found in the weak- 
coupling BCS limit [35]. 

To derive the effective induced interactions among the 
spins-t, mediated by the spins-J,, wc derive self-consistent 
equations for the corresponding Green's functions for the 
two species of fermions. As discussed above, we assume 
the spins-J, to be in an unpaired Fermi-liquid state, while 
the spins-t may possess pairing correlations. For the 
latter, accounting for such pairing correlations is best 
done by organizing the spin-f Green's functions using 
the Nambu notation: 



*t« - ( ^i|:J ; • 



(3) 



We can then define the imaginary time ordered Nambu 
matrix Green's function for spins-f 

G^(r,r)--(r$^(r,r)<i>j(0,0))^, (4) 

with matrix elements 

The normal Green's functions arc given by 

G^(r,T)--(rvl/^(r,r)*|(0,0))^, (6) 

and the anomalous ones are 



F^(r,T) = -(T*^(r,T)*^(0,0))^, 
Ff(r,T) = -(TvI/t(r,r)M'|(0,0))^. 



(7) 



For notational convenience, we further introduce the 
four- vector k = (k, ia;„). In Fourier-Matsubara space 
the spin-t Green's function satisfies a matrix analog of 
the Dyson equation 



Gt(fc) = Go,t(fc) + Go,t(fc)St(fc)Gt(fc), 



(8) 



where 



Go,t(fc) 



Go.t(fc) 

-Go,t(-fc) 







iwii+skt 



with ^ko- = Ck - A*(T! and 



St(fc) 



Et(fc) At(fc) 
A;(fc) -I]t(-fc) 



(9) 



(10) 



is the matrix self-energy. The appearance of nonzero off- 
diagonal terms in the self-energy indicates the presence of 
a supcrfluid state, with order parameter or gap function 
A|(fc). The self energy satisfies 

St(fc) = XCT, Y. ^i(9) + ^'^- E Gt{q)Gi{q')x 

q q.q' 

X r(g, q', q + q'-k, k)Gi{q + q' - k), (11) 

where (X^ is a Pauli matrix, T is the reducible two-particle 
matrix vertex function (apart from energy and momen- 
tum conserving delta functions), and the summation is 
E, = (/3V)-^ EqEi.., with /3-1 = k^T and V the sys- 
tem volume (which, henceforth, we will set to unity). 
Diagrammatically, Eq. (11) is shown in Fig. 2. 

The spin-down Green function Gi (r, r) satisfies a sim- 
ilar set of equations (but is assumed to be unpaired) ; in 
terms of the spin- J, fermion operators it is defined as 

G^(r,r) = -(TvI'^(r,r)vI/t(o,0))«. (12) 

Since the Greens's function depends on "^"[{k) via the 
Dyson equation [Eq. (8)], this, along with Eq. (11) 
amount to self-consistent equations for the self-energy. 
For pairing to be stable, we must find a solution to these 
equations possessing a nonzero off-diagonal component 
A-|-(fc). Formally, the only assumption we have made 
thus far is that the spins-4. are unpaired (i.e., they pos- 
sess no off-diagonal component to their self energy); in 
practice to proceed we must make a physically-motivated 
approximation for the Bethe-Salpeter equation satisfied 

byr. 



III. LEADING-ORDER PERTURBATION 
THEORY 

Having set up the general formalism for computing the 
self-consistent pairing amplitude for intraspccics pairing 




St(fc) 



FIG. 2: The formally exact diagrammatic matrix self-energy, 
Eq. (11). The spin-t Green's function and vertex function 
r are Nambu matrices, while the spin-^ are normal scalar 
functions. The interaction lines are Xer^, where er^ is the 
third Pauli matrix. 



correlations among the spins-'f, in the present section we 
take the weak-coupling pcrturbative approximation A -^ 
0— in which these equations simplify. This will allow 
us to easily identify the effective interaction among the 
spins-t; as we show below, a similar structure will hold 
when we sum diagrams to all orders in A (to access the 
unitary regime) . Our results in this section are consistent 
with those of Ref. 35. 

In the weakly interacting BGS limit we can expand per- 
turbatively in A, or equivalently in Os = r?iA/47r. In fact, 
to obtain a nonzero result for the intraspecies pairing we 
must keep terms of order A^. The reason for this is that, 
to 0{\), the self energy is simply given by the first term 
in Eq. (11), the Hartree contribution, H^ik) = Acr^nj^, 
where n^ is the density of spins-4-. At this level of ap- 
proximation, the Nambu self-energy is diagonal and no 
pairing is induced. 

To obtain results valid to quadratic order, 0{\^), it is 
sufficient to approximate the vertex function by the bare 
interaction. 



-Act, 



(13) 



where the presence of the Pauli matrix is due to the fact 
that we need to write the interaction in Nambu space. 
This gives for the self energy: 

St(fc) = Acr^?i^ - A^<T^ y^ G-x{q)(TzGi{q)G^{q + q - k). 

q.q' 

(14) 
Here and in the next section (when we sum diagrams to 
all orders in A), we'll neglect the diagonal component of 
the self-energy (simply assuming it amounts to a chemical 
potential shift). Within this simplifying approximation, 
the Green's function for the spins-f can be written as 



Gt-(k, ia;„) 



-w„ + ^kt At-(k,iw„) 

{iuJnf -El\ A:f (k, iuJn) ioJn - ^kf 



(15) 
where £;k(ia;„) = W^^^ + A|(k,ju;„)A^(k, ia;„). Gom- 

paring Eqs. (5) and (15) one can read off the relationship 
between the anomalous propagator and gap function 



F-t.(k,zw„) 



At(k,ia;„) 



w„ 



K 



(16) 




FIG. 3: The induced interaction among the spins-t, mediated 
by the spins-J,, in the weak-couphng BCS hmit. The internal 
loop corresponds to density fluctuations of the minority spin. 
The external legs are only present to label the incoming and 
outgoing energy and momentum. 




FIG. 4: (Color online) The p-wave projection of the Lindhard 
Function, see Eq. (20), as a function of density imbalance z = 
fcp^/fcpj,. The peak near z ~2 occurs as a result as the turning 
off' of particle-hole excitations with a transferred momentum 
larger than the diameter of the spin-4. Fermi surface, while at 
smaller wave vectors the Lindhard function is approximately 
ffat, leading to a vanishingly small p-wave projection. 



This approximation also leads to a simple form for the 
upper-right off-diagonal component of Eq. (14): 

A-,.(k,iw„) = -^2^ V;ff(k,q,ia;„,ii/„)^^2-rT^' ^^^> 



q-^t^n 



n ' q 



Feff(fc,g) = A2^G4(fc-g + p)G^(p), (18) 



where, in the second line, we used a short hand notation 
T4ff(fc, g) = V^fj(k, q, icj„, ii'n). Thus, the similarity of 
the resulting expression to a standard gap equation for 
pairing has allowed us to identify an effective interaction 
T4ff(fc, fc') that is plotted diagramatically in Fig. III. 

We proceed by approximating the full spin- J, 
Green's function in (18) by noninteracting ones, i.e., 
G'4,(k, iwn) -> Go,;(k, iw„) = (iw„ - ^k4,)~\ neglect- 
ing the frequency dependence of T4ff(A:,fc') (setting the 
Matsubara frequencies to zero), and assuming that 
A-t-(k, iw„) = A-f(k), i.e., it is indepenent of frequency. 

Evaluating the remaining Matsubara sum in Eq. (17), 



we arrive at the following gap equation that is of the 
standard form: 



At(k) 



,,-^ tanh (/?£:„ /2) 
5^yeff(k,q) A^^i^A^(q), 

q q 



(19) 



where now the induced interactions (mediated by den- 
sity fluctuations of the spins-J,) are proportional to the 
Lindhard function: 



Feff(k,q) = - 



4:71 a. 



7V4(eFjL(|k-q|/2fcFi))(20) 



4x 



In 



1 



(21) 



where we have replaced A — >■ Anas/m as noted above 
(valid in the weak-coupling BCS regime). Here, fcpo- is 

the Fermi wavevector for species a, satisfying rio 
and 



3ir2 



7V^(eFa) = mfcF^/(27r2), 



(22) 



is the density of states at the spin-tr Fermi surface. 

We now decompose the gap potential and effective in- 
teraction into angular momentum channels using 



At(k) 



E(2^ 

oo 

E(2^ 



l)v'Pi{k-k'), 



1)A' Pi(z-k), 



(23) 
(24) 



1=0 



where Pi{x) are the Legendre polynomials. The transi- 
tion temperature for each angular momentum channel is 
then determined by the solution to 



, _ 'r- , tanh (/3eqt/2) ^; 

^kt - - 2^ ^k,q ^^— ^9t' 



(25) 



where we have assumed a continuous transition where the 
gap potential vanishes; in this limit we replace iJq — s- S,q\. 
To proceed, we note that we are interested in the onset 
of p-wave {I — 1) pairing at the spin-f Fermi surface. 
Thus, we set k -^ k-p^ in Eq. (25), and henceforth choose 
I — 1. The summation over q is then dominated by 
the region where q ~ k^^- Using that the function v\ 
is only nonzero for k and q within fcpj, of each other, 
and converting the sum to an integral (introducing the 
density of states), we obtain 



A|(eF^) w -4pi.,fepT^T(eFt) 



tanh (^e/2) . , , 
de ^f^^Ai(6), 



(26) 
for the pairing gap near the assumed continuous transi- 
tion. The corresponding transition temperature is: 



kv.Tr 



2e 



F4,e' 



■ exp 



1 



«fcFt,fcFt^t(eFt) 



(27) 



where 7 = 0.577... is the Eulcr-Mascheroni constant, 
giving our result for the transition temperature expressed 
in terms of the £ = 1 projection of the induced interac- 
tions in the weak-couphng regime a^ — > 0~. 

The final step in the perturbative analysis is to obtain 
the p-wave component of the perturbative effective inter- 
action, Eq. (20). Then, wc have for the dimensionlcss 
effective induced interaction [appearing in the argument 
of the exponential function of Eq. (27)] [35]: 



where z = fept/^Fj. a-nd 



-iikpiasf 



L\z) 



5z^ 



15z4 



■In 1 



2^ + 5 

30z 



In 



zL\z), 



1 -z 



l + z 



is the p-wave projection of the Lindhard function 



7r 

L^{z)= d9 coaeL{\k-k'\z/2), 



(28) 



z^ + 2 

T5^' 
(29) 



(30) 



T(fc) = 



A + 



FIG. 5: The T-matrix, which describes the ladder sum of 
two particles interacting through the bare potential infinitely 
many times. Here, the solid lines are Greens functions and 
the wavy line is the interaction A. 



the strongly interacting regime, near unitarity, one must 
sum an infinite number of diagrams, to all orders of the 
bare interaction A, contributing to the Nambu self-energy 
Eq. (11). We begin with the conventional T-matrix ap- 
proximation, which, when extended to a system with 
pairing correlations, includes ladder and crossed-ladder 
diagrams. 



A. Ladder plus crossed ladder approximation 



with k ■ k' — cos{9). The perturbative formula for Tc, 
given by inserting Eq. (28) into Eq. (27) thus yields a re- 
sult that vanishes exponentially as kBTc exp[— c/(fcFas)^] 
with c > 0; if this estimate is correct, then such p-wave 
pairing is probably not experimentally observable in the 
weak-coupling perturbative regime. 

For fixed kpiUs , the density imbalance dependence en- 
ters through zL^{z), that we plot in Fig. 4. We see that 
this quantity shows a maximum value near z « 2, lead- 
ing to a peak in the p-wave transition temperature near 
a polarization of P w 0.77; we find a similar peak in Tc 
in our unitary-regime results to follow. 

In the next section we proceed to derive a formula for 
Tc which goes beyond the weak-coupling regime. Our 
result is of the same form as Eq. (27) but with a more 
complicated expression for the effective induced interac- 
tions vl I, , including contributions from all orders of 
A. Our final result for Tc in fact reduces to the perturba- 
tive result of this section when we take the weak-coupling 
limit as — > 0—; however, in the regime where they agree 
(far to the right of the displayed area of Fig. 8), Tc is 
many orders of magnitude smaller than the maximum 
transition tcinperature shown in Fig. 1. 



IV. BEYOND LEADING ORDER 

In the preceding section, wc showed that, within 
leading-order perturbation theory, there is a transition 
to a p-wave supcrfluid of the majority species of a 
population-imbalanced Fermi gas, although the predicted 
perturbative temperature is vanishingly small. In the 
present section, we show how additional classes of dia- 
grams, occurring near the unitary regime, can lead to an 
enhanced Tc- To investigate the induced interaction in 



Recent work has found that the nonsuperfluid Fermi 
liquid phase of imbalanced Fermi gases is well described 
by the so-called T-matrix approximation for the self- 
energy or more specifically for the vertex [26, 27, 32], 
depicted diagrammatically in Fig. 5. In the case of un- 
paired spins-f, this approximation amounts to summing 
the repeated interaction of a spin-f particle with a single 
spin- 4, particle-hole bubble, as shown diagramatically in 
the self energy shown in Fig. 6(a). 

Our aim is to generalize the T-matrix approximation 
to imbalanced Fermi gases by including the possibility of 
pairing for the spins-f. The simplest way to do this is to 
keep the same series of diagrams, but replace the normal 
Green's function G'^(fc) with the corresponding Nambu 
Green's function, i.e., considering the same ladder series 
but with Gf{k) — > G-|-(fc). We refer to this Nambu T- 
matrix with the bold symbol T(fc). Formally it is the 
same as Fig. 5 but with the spin-t Green's functions pos- 
sessing Nambu structure (while the spin- J, Green's func- 
tion is still a normal Green's function); additionally the 
coupling (wavy line) is given by Acr^ . The corresponding 
self energy is given in Fig. 6(a). 

However, this ladder series is not sufficient for our pur- 
poses, as can be seen by noting that it docs not even 
reproduce the standard T-matrix approximation in the 
limit A-|-(fc) — > 0. This is because, in the Nambu nota- 
tion, the lower right element of Eq. (5) can be regarded 
as a Green's function with a line possessing opposite mo- 
mentum from the upper left clement, i.e., its momentum 
is in the opposite direction. This means that, in addition 
to ladder diagrams, we must include maximally-crossed 
diagrams shown in Fig. 6(d). 

The conclusion of the preceding remarks is that the 
set of diagrams needed to minimally generalize the usual 
T-matrix approximation to a system with pairing among 
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FIG. 6: The diagrammatic contributions of the spin-t Nambu 
self-energy at various levels of approximations. Where T, de- 
fined by Eq. (32), corresponds to the Nambu matrix gener- 
alization of the T-matrix, shown in Fig. 5, where all spin- 
t Green's functions are in Nambu space while the spins-^ 
are scalar functions and interaction lines correspond to Acr^. 
Here Tc, given by Eq. (33), is the Nambu generalized, sum of 
the so-called maximally crossed diagrams, shown in panel d, 
T = HaT and Tc = -Ug^Tc, with Ha and n^ defined by 
Eq. (34) and Eq.(35) respectively. 



the spins-f is given by Fig. 6(b). The first term contains 
the ladder series and the second term contains the crossed 
diagrams Fig. 6(d). This gives: 

E^(fc) = 5] T(fc + q)G^iq) + J2 Tc(g - k)G^iq), (31) 



where 

T(fc) = ^4lA-i+n^(fc)]"\ (32) 

T,{k) = AV,n|(fc)[lA-i+nB(fc)]"' (33) 

Here, the inverse is understood as a matrix inverse and 
the bubbles IlA{k) and Usik) are given by 

n^(fc) = J2G^{q)G^ik-q)cT,, (34) 

Usik) = J2G^iq)G^ik + q)a,. (35) 



Having derived the natural generalization of the lad- 
der approximation to incorporate the possibility of in- 
traspccics pairing among the spins-f, we now show that, 
in the unitary regime, it cannot possess any ofF-diagonal 
pairing correlations but merely reproduces the polaron 



Fermi liquid, implying that we must go beyond this level 
of approximation. To do this, we make the replacement 
Eq. (2) and take the limit A — > c« while holding Ug 
fixed (a procedure that, above, we argued to be valid 
away from the weak-coupling BCS limit). The Nambu 
T-matrix is thus 



T(fc) = (T, 



47ra, 



n.(^)-iEi 



2ek 



(36) 



where, when we take the limit A — > oo, it is clear that the 
last term yields a divergence proportional to the identity 
matrix. To be nonzero, a similar divergence must appear 
in n/i(fc); however, because of the Nambu structure of 
G-|-(g), only the upper left component of 11^ can possess 
such a divergence. Thus, all elements of T(fc) except the 
upper-left will vanish, and we obtain: 



lim T(fc) 



_ f T(fc) 




(37) 



where T(fc) is the usual T-matrix, satisfying 
Anas 



[T(k, iuj,i)] 



lUJr, — IVr 



(38) 
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To be clear, the quantity G'^(q, i;/„) appearing in this 
formula is the upper-left component of the full matrix 
Green's function G^(q, w„). A similar simplification oc- 
curs in Tc(fc), which, when we make the same replace- 
ment and take the hmit of A — >■ oo, yields 



lim Tc(fc) 



giving for the self-energy 
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st(fc) = Y. 



J{q + k) 





-T(g- 



k) 



(39) 



G,(g), (40) 



which, crucially is diagonal in Nambu space and, thus, 
possesses no pairing correlations. In fact, Eq. (40) is ex- 
actly the self-energy within the usual "polaron" Fermi 
liquid description of strongly imbalanced Fermi gases, 
equivalent to the Chevy variational wavefunction [30] as 
shown in Ref. 31. 

Thus, although we generalized the ladder approxima- 
tion in the simplest possible way to include pairing corre- 
lations (summing ladder and crossed-ladder diagrams but 
with a Nambu spin-t Green's function), wc have found 
that, within this approximation, such pairing correlations 
are not stable. The results of this section, however, tell 
us two things: Firstly, we must consider sets of diagrams 
that go beyond this approximation. Secondly, to have a 
result that is finite when we make the replacement Eq. (2) 
and take A — >■ oo, we must consider diagrams with ladder 
or crossed-ladder type subdiagrams. 



B. Combined ladder and crossed-ladder diagrams 

As we have seen, the ladder and crossed-ladder (or 
maximally-crossed) sets of diagrams contain geometric 
sums that are nonzero when we exchange the bare cou- 
pling A for the scattering length as using Eq. (2) and 
take the limit A — > oo. Since the quantities T(fc) and 
Tc(A:) arc diagonal in this limit, each possessing only 
one nonzero element, the corresponding contributions 
to the self-energy in Fig. 6(b) are also diagonal. How- 
ever, there exist additional sets of diagrams containing 
T(fc) and Tc(fc) as subdiagrams that are not diagonal 
in this limit; the simplest such self-energy diagrams arc 
shown as the final two diagrams in Fig. 6(c). Our inclu- 



J 



sion of these diagrams may alternatively be understood 
on physical grounds as due to the fact that paired su- 
pcrfludidity mixes particles and holes: roughly speak- 
ing, T describes particle-particle scattering and Tc de- 
scribes hole-hole scattering. Including pairing requires 
mixing these, suggesting the incorporation of these di- 
agrams. Thus, Fig. 6(c) represents the full set of self- 
energy diagrams that we consider here: The ladder and 
crossed-ladder terms that yield the polaron Fermi liquid 
(as shown above) and the combined ladder and crossed- 
ladder diagrams that, as we now show, capture the insta- 
bility of the polaron Fermi liquid to p-wave pairing. This 
self energy is: 



S^(fc) = 5^ T(fc + q)G^iq) + J2 Tc(g - k)G^iq) 

- AV, Y. ^A{k + q) [lA-1 + n^(fc + q)] ~^G^{q)(TMB{q - q) [lA"! + Ueiq - q)] ~^G^{q' - q + k)G^{q') 

- AV, Y1 ns(9' - k) [lA-1 + UA{q' - k)] ^'Gt(g)<T,n^(g + q') [t\-^ + IlA{q + q')] ~^G^{q + q' - k)G^{q'). 

9,9' 

(41) 



Again exchanging A for a^ using Eq. (2) and letting A 

oo, we get 



Acr,n^(fc)[lA-i+nA(fc)] ' 



J{k) 




and 
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S^(fc) = y I "^(« + ^)^i(«) V^g{k,q)F^{q) 



A(T,nB(fc)[iA-i + nB(fc)] ' 

so that Eq. (41) reduces to 

T(g + fc)( 

V^^ik,q)F;iq) -T(q-k)G^{q) 

^ . (44) 

where the effective pairing interaction V^g{k, k') is shown 
diagrammatically in Fig. 7 and is explicitly given by 

V;ff(fc, fc') =J2T{q + k)T{q - k')G^{k -k' + q)G^{q). 

(45) 
The effective interaction Eq. (45) is a direct generaliza- 
tion of the induced interaction of the weak BCS limit 
(18), except now with a energy and momentum depen- 
dent coupling, i.e., the bare interaction A in Eq. (18) has 
been replaced by the T-matrix in Eq. (45). 

With this effective interaction, the gap equation is 
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FIG. 7: The effective induced interaction between spins-t, in- 
cluding contributions to all orders in perturbation theory, that 
is the many-body generalization of the leading order effective 
interaction shown in Fig. Ill 



where, as in the weak-coupling case, in principle A-|-(fc) is 
frequency-dependent, although we shall again assume it 
to be static. Unlike the weak-coupling regime, however, 
now the effective interaction also depends on the pairing 
gap via the relationship, Eq. (38), of the T matrix to the 
spin-f Green's function. This difficult self-consistency 



problem simplifies near the assumed continuous phase 
transition, where A-|.(k) vanishes, a regime we now focus 
on. In this regime, we can simply approximate the T 
matrices by their form in the imbalanced normal phase. 
In principle all of the Green's functions appearing in 
the effective interaction Eq. (45) are the exact interact- 
ing Green's functions. To proceed, we make some sim- 
plifying approximations, firstly by replacing these with 
their noninteracting expressions (equivalent to assuming 
the imbalanced Fermi liquid phase is noninteracting) and 
neglecting the external frequency dependence of Eq. (45) 
(i.e., setting cj/; = w[. ^> 0). This gives 

Kff(k,k') = 2r^y^T(q + k,zr!)T(q-k',*f)) 

1 1 

(47) 
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The sum over Matsubara frequencies can be evaluated 
using the standard trick [44] that requires the location of 
the poles of the sumniand in the complex plane. This is 
simplified by our knowledge that poles of the normal- 
state T-matrix on the real axis occur at the onset of 
s-wave paired superfluidity. Since we are studying the 
strongly imbalanced limit, we can assume no such poles 
contribute, and proceed by keeping only the poles from 
the Green's functions in the second line of Eq. (47). This 
leads to 

Kff(k,k')-2r^y'^T(q-Hk,eq4)T(q-k',eqi) 



"F(Cqi) 



(48) 



our final result for the static induced interaction for the 
majority spins-t, mediated by the spins-J,, in a strongly 
interacting imbalanced Fermi gas. 

Our main results come from numerically evaluating 
Eq. (48). However, it is useful to first make a bit more 
approximate analytic progress by invoking an approxi- 
mation that is valid in the large-imbalance limit. Thus, 
to evaluate the remaining integral over q, we begin by 
noting that the Fermi function restricts, at low T, the 
momentum to q < fcp^, inside the spin- J, Fermi surface. 
But since we're interested in the regime where k and k' 
are on the spin-^ Fermi surface, in the strongly imbal- 
anced limit k-p-f ^ k^i the momentum argument of the 
T-matrices in Eq. (48) is approximately given simply by 
/cp^. In this strongly imbalanced limit the first T-matrix 
in Eq. (48) simplifies to 



T(kF,eqi)=^( 



mk 



Ft 



47ra, 



47r2 



(49) 



with a similar expression holding for the second T-matrix 
in Eq. (48). Within this approximation, valid at large 
population imbalance P — > 1, the T-matrices are thus 
independent of momenta. As can be seen in the origi- 
nal expression Eq. (44), when this occurs T4fF is simply 



P 



0.8 




0.6 




0.4 
0.2 
n n 


"S. 

1!^^^^^^^ .1.1. 



-0.2 0.0 0.2 0.4 0.6 O.i 

-{kpas)'^ 



1.0 



FIG. 8: (Color online) The p-wave channel of the effective 
interaction between majority spins, vl j. , multiplied by 

the Fermi-energy density of states Nf{e-p^), as a function of 
the density imbalance P and s-wave scattering length as, is 
shown. Here, kp — {k-pf + fcp^)/2. At large P, above the 
dashed line, it is attractive leading to a p-wave superfluid 
at temperatures below T^. The solid white line labels the 
location of a line of FFLO quantum critical points, defined 
by T(Qfflo,0) — >■ 00, at zero temperature. This coincides 
with the location where vl j. is large. 



proportional to the bare bubble occuring in the weakly 
interacting regime (see Eq. (18)). Evaluating the remain- 
ing integrals and making the projection to the p-wave 
channel (as we did in the weakly interacting regime), we 
obtain 



Kg Jc 



2e'> 



-ep^exp 



2 ln{z) 






(50) 



with z — kp^/k-p,, a result that we emphasize is only 
valid in the asymptotic strongly imbalanced regime. 

Away from this P — > 1 limit, we must perform a nu- 
merical evaluation of the integral in Eq. (48) and the 
p-wave projection; this yields a dimensionless effective 
interaction v^ j, iV-|-(eF-|-) that we plot in Fig. 8 and in 
the top panel of Fig. 10. For this calculation we assumed 
the r — >■ limit; since our resulting transition temper- 
ature is still small compared to the spin-f Fermi energy 
this should be an accurate approximation. These figures 
show that we find attractive interactions in the p-wave 
channel over a wide range of the interactions and the 
population imbalance. 

The corresponding transition temperature, plotted in 
the bottom panel of Fig. 10 as well as in the phase di- 
agram Fig. 1, is given by the same formula as in the 
weak-coupling limit, i.e.. 



knTr 



2e 
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where to arrive at this formula we used the fact that, as 
in the weak-couphng hmit, v^ j., (the p-wave projection 
of Eq. (48)) is only nonzero for k and k' close to each 
other, within a window of approximately ±ep, . We have 
verified this numerically; a plot displaying the typical 
behavior of v], j,, is shown in Fig. 9. 

It is interesting to compare the value of Tc obtained via 
a direct numerical evaluation of Eq. (48) to our approx- 
imate asymptotic formula Eq. (50). In Fig. 11 we plot 
these as a function of population imbalance for three dif- 
ferent values of the scattering length. These curves show 
that, while Eq. (50) is accurate asymptotically close to 
P — 1 , it misses the peak and drop in T^ that generically 
occurs with decreasing P. 

What is the origin of this peak in the predicted T^ 
occuring for very large imbalance? Certainly, one ex- 
pects an increase in T^ as P decreases from unity, arising 
from the increase in the density of spins-J, which provide 
the induced interactions. More precisely, the effective 
interactions among the spins-f are mediated by particle- 
hole excitations of the spins-^. However, particle-hole 
excitations with a wavevector larger than the diameter 
of the spin-down Fermi surface (2fcF4,) are energetically 
suppressed, implying that the density response function 
is strongly varying for k ~ 2A:f^. This strong variation 
as a function of momentum leads to a large p-wave pro- 
jection of the induced interaction when k ~ ky^ (since, 
for p-wave pairing, we require an attractive interaction 
that strongly varies around the spin-^ Fermi surface). If 
we then assume that the maximum T^ will occur when 
/cpt = Sfcpj^, we are led to the prediction that Tc will peak 
near 
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close to the value of the peak position shown in Fig. 10. 

We note that Eq. (52) only approximately locates the 
peak position; indeed. Fig. 10 shows some variation of 
the peak position as a function of interactions. (We have 
not pushed this computation deep into the BEC regime 
where we know the magnetic superfluid ground state in- 
tervenes [12, 15].) Indeed, the preceding argument is 
strictly true in the weak-coupling limit where the effec- 
tive interaction is given by the Lindhard function L{x) 
that has a singularity near a; = 1 that leads to a sim- 
ilar peak in T^ in the weak coupling limit. (Although, 
as mentioned above, the weak-coupling Tc is orders of 
magnitude smaller than the values plotted here.) In the 
strong-coupling limit, we expect this argument to still 
approximately hold since there the Luttinger theorem 
ensures that the presence of spin-'t" and spin- J, Fermi sur- 
faces at the same volume as in the weakly interacting 
limit (until the broken symmetry phase appears) [45, 46]. 

We also note an additional possible reason for the oc- 
curence of a peak in Tc at large imbalance: A proximate 
FFLO instability. Thus, one expects an instability to- 
wards FFLO pairing for imbalanced Fermi gases, occur- 
ing when there is a divergence of the retarded T matrix at 
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FIG. 9: (Color online) The momentum dependence of the 
p-wave channel v], f., of the effective interaction, Eq. (45), as 
defined by Eq. (23), for zero frequency and at unitarity for a 
polarization of P = 0.85. As one can see the interaction is 
nonzero only within a window of approximately fc = fc' ± fcp^ . 



zero frequency (but nonzero wavevector Qfflo)- In the 
vicinity of such a phase transition, the T matrix would, 
correspondingly, possess a large magnitude that would 
enhance the effective induced interactions. To test this, 
in Fig. 8, we plot, as a white line, the P at which such a 
zero temperature FFLO instability would first occur; as 
seen in this plot it closely correlates to the regime where 
the p-wave Tc is largest. Our calcultion of this FFLO 
transition of course neglected the possibility of p-wave 
pairing among the spins-f. Thus, if the p-wave transition 
occurs first, with decreasing temperature, it would likely 
move the location of the FFLO phase boundary. Further 
detailed analysis will be required to sort out these various 
competing instabilities. 

We conclude this section by briefly justifying the use of 
the on-shell approximation, in which we assume that we 
only need the p-wave induced interactions for momenta 
on the spin-f Fermi surface. As is known from the theory 
of superconductivity, such an approximation is typically 
only valid for weakly interacting quasiparticles and, near 
unitarity, particles of opposite spin are strongly interact- 
ing. Despite this, in highly imbalanced systems the life- 
time of the quasiparticles remains extremely long; this 
implies the quasiparticle-quasiparticle interaction is still 
quite weak. This can be seen directly from the experi- 
mental and theoretical results of Ref. [28] . In Ref. [28] the 
Fermi liquid properties of systems with an imbalance as 
low as P « 0.70 were explained well within a theoretical 
framework consisting of only a single spin- J, quasiparti- 
cle, which indicates induced interactions of like spins, is 
subdominant to such things as the rcnormalization of the 
chemical potential and effective mass, at least for T >Tc. 
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FIG. 10: (Color online) The top panel shows the on-shell £ = 1 
channel of the effective interaction (45) for spin-f fermions, 
times the density of states, as a function of polarization, at 
unitary and into the BEC side. The bottom panel shows 
the corresponding p-wave transition temperature, according 
to Eq. (27), with ept/ep.^ = [(1 + P)/(l - P)]'''^ 
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FIG. 11: (Color online) A comparison of the transition tem- 
perature obtained from Eq. (50), vaUd for asymptotically 
large polarizations, and a full numerical treatment of the in- 
duced interaction Eq. (45), using Eq. (51), with ept/eFj. ~ 

[(l + p)/(l_p)]2/3. 



Thus, wG argue that our analysis is valid if the strongly 
imbalanced nonsuperfluid phase is truly a Fermi liquid 
for T>T,. 



V. LOW-TEMPERATURE PAIRING GAP 

In the present section our goal is to estimate the mag- 
nitude of the p-wave pairing gap at T — > 0. One subtlety 
is that our method for calculating the induced attrac- 
tion among the spins-f is only valid for T > Tc, since 
we neglect the pairing among the spins-f in computing 
the T matrix. However, if we assume that the induced 
pairing A-|-(k) is small in magnitude compared to epf (as 
one might expect), then we can assume A-|-(k) has only 
a small effect on the T matrix, and proceed to neglect it 
when estimating A-|-(k). 

The relevant equation for the low-temperature gap is 
Eq. (19), but with ycff(k, q) given by Eq. (48). Focusing 
on the p-wave channel, we have, for T — > 0, 
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Following the work of Anderson and Morel, we expect 
the dominant p-wave instability to be of the Px + ipy 
form [47]. We thus write 



A^(k) =: A(fc2)(fca; +iky 



(54) 



Converting the sum in Eq. (53) to an integral (recall we 
set the system volume to unity), we obtain 
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arises from evaluating the angular integration. We now 
use the fact that vl. is sharply peaked near k = q, and 
assume the rest of the integrand of Eq. (55) is smooth 
there. This yields the approximate formula 
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Restricting attention to the vicinity of the Fermi surface 
by setting k — )■ fcp-|-, and using the asymptotic value of 
F[x\ -^ |x for X 



A(fc2 



0, we have 
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(59) 



We can further simplify this formula by noting that, since 
v], is peaked for k = q with a width of approximately 
ifcpj,, we have Vk^^ ~ 2kF^vl^^j.^^. Then, normalizing 

and using Eq. (22), 

N^{eFty{kF^,kFt)- (60) 



A(fcp) to the Fermi energy and using Eq. (22), we have 
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The product of the last two factors of Eq. (60) is pre- 
cisely what is plotted in Fig. 10. At unitarity, aj^ = 0, 
this factor reaches ~ 0.4 at P ~ 0.85, or ^ ~ 0.43. 
Plugging these values into Eq. (60) gives the estimate 
^ pa 0.4. Although this estimate is rather large, it 
is important to keep in mind that this is the maximum 
pairing gap; the full gap function will exhibit nodes ac- 
cording to Eq. (54). An important issue for future work 
is to find a more accurate estimate for the low-T pairing 
gap for the majority species of an imbalanced Fermi gas. 



VI. EXPERIMENTAL DETECTION 

We argue that, although small (a few percent of the 
Fermi energy), the transition temperature is within range 
of current experimental capabilities [55]. However, the 
experimental detection of such a state may still be chal- 
lenging. One striking way to identify the presence of 
superfluidity is via the presence of vortices in a rapidly- 
rotating cloud, as done in Ref. 51 to detect s-wave pairing 
correlations. However, in Ref. 51 the detection of BCS 
pairs was accomplished by ramping the magnetic field 
onto the BEC side of the resonance, changing them into 
molecular pairs. For the case of the p-wave Cooper pairs 
predicted here, it is not clear that such a ramp is possible. 

Owing to this difficulty, here we focus on two other 
possible ways to detect px + ipy Cooper pairing among 
the majority species of an imbalanced Fermi gas: Radio- 
Frequency (RF) spectroscopy (in which pairing correla- 
tions are detected by a shift in the rate at which atoms in 
a particular state are transferred to a third unoccupied 
level [48]) and noise correlations (studied theoretically in 
Ref. 49 and experimentally in Ref. 50). 



A. RF spectroscopy of the majority spin p-virave 
superfluid 

To use RF spectroscopy to probe Px + ipy pairing of the 
spins-t, one would measure the number of spin-f atoms 
transferred to a known unoccupied level, as a function of 
applied RF frequency [53]. Neglecting final state effects, 
the total transferred at a given frequency is given by 

I„{uj) oc^A^(k,4^ -^)nF{£.ua ~^), (61) 

k 

where the spectral function A^^ (k) is related to the imag- 
inary part of the retarded Green's function G^^{k) by 
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A,(fc) = --ImG;(fc), 



(62) 



and np(w) is the Fermi distribution function. It is con- 
venient to write the superfluid Green's function in the 
mean-field spectral representation: 
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FIG. 12: (Color online) The normalized zero temperature 
RF line-shape (61), I{i^) — I^(uj)/ j?° du) I^iui), of a p-wave 
superfluid, with a. px + ipy ground state. The magnitude 
of the gap is approximated by the transition temperature, 
Ao ~ 0.03ep.|. and the chemical potential p.^ ~ ep^, for 
P ~ 0.85 at unitarity. In general for a. px + ipy ground state 

the peak occurs near a; ~ i/M? + ^o ~ /^ti ^^'^ as a; — )• oo, 

I{cj) -^ i^-'^l'^, while for cj ^ 0, /(lj) ^ (o; + 2ii^fl'\ 



where the so-called coherence factors are 
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GUKlu) 



UJ - Ei^ + ii] uj + Ey + irj 



(63) 



(64a) 



The spectral function is 



vl5{uj + E^). (65) 



As we have discussed, the p-wave ground state is expected 
to have px +ipy symmetry, given by Eq. (54) [40, 47, 54]. 
A full calculation of the RF lineshape requires A{k'^) for 
all k, a difficult self-consistency problem that is beyond 
the scope of this work. We proceed by simply assuming 
A(fc2) = Aq, i.e., a constant value. At zero tempera- 
ture and using Eqs. (64a), (65), and (54) the momentum 
integrals in Eq. (61) can be done exactly, although the 
result is too unwieldy to present here. Figure 12 shows 
the resulting RF line shape at unitarity. In this plot we 
chose Aq ~ fcel^ci assuming the magnitude of the pairing 
gap reflects the transition temperature (with the latter 
given by Tc ~ .03eF^/fcB! a typical maximum value in 
Fig. 10). Such an estimate is more conservative than the 
rather large estimate found in Sec. V. As can be seen from 
Fig. 12, unlike the s-wave state, which has a hard gap for 
u! < Aq/cf^, the RF line shape of this p-wave state re- 
mains non-zero for all uj, due to the nodes in the pairing 
gap. In principle one could use this to detect the p-wave 
phase of imbalanced Fermi gases. One issue, however, is 
the small magnitude of the peak position in the energy in 
Fig. 12. which occurs at a scale of order Aq/ep-^ (rather 
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than at Ag). To see whether p-wavc pairing is truly ob- 
servable, we need a good estimate of the low-temperature 
pairing gap in this phase. 



B. Density correlations 

The use of spatial correlations in the density of the 
free expanded gas, or density-density correlations, as a 
tool to probe the order of the superfluid state of ultra- 
cold atomic systems was put forth in Ref. 49. As free 
expansion of the density probes the momentum distribu- 
tion of the trapped system, spatial correlations of this 
expansion probes correlations in momentum space. It is 
these momentum correlations of the superfluid state that 
are a direct consequence of the Cooper pairing. 

Theoretically the quantity of interest is the equal time 
density-density correlation function of the spins-f 

D{v, r', t) = Trpt,,p 5n^{v, t)5h^{v' , t), (66) 

where Sh^{r) = n-|-(r) — {h^{r)) and Pt^i^p is the equilib- 
rium density matrix of the trapped system, i.e., super- 
fluid, while the time evolution of the operators is given 
by free un-trapped Hamiltonian. As the spins- J, arc as- 
sumed unpaired, their momentum correlations arc essen- 
tially featureless and won't be included. Neglecting the 
inhomogeneity caused by the trapping, we take the den- 
sity matrix to be given by the zero temperature BCS 
state. 



I* 



BCS 



)-n 



^k ^l^;-]- '^ _ kl" 



I vac), 



(67) 



with the coherence factors given by Eq. (64a) above. 
Within this approximation and assuming the density is 
measured on space and times scales much larger than fcjL, 
and ep|, we find 

D{r, r',t)^ (£^) ' J2 0* (rm/t)0q(r'm/t) 
q.q' 



X (?iq,(r'm/i)(/)*,(rm/t)|uqp|z;q| 



(68) 



where 0q(k) is the Fourier transform of the single-particle 
wave function of the trapped system with quantum index 
q. Typically these will be very sharply peaked near k and 
k', thus 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



-1.5 -1.0 -U.5 0.0 0.5 1.0 1.5 
(b) 



FIG. 13: The column integrated density-density correlations 
(70) (weight of delta function only) of a free expanded Px+iPy 
superfluid with maximum gap Aq = 0.03ep.|. and symmetry 
axis z is shown: (a) by integrating along the z-axis and (b) 
integrating along the j/-axis. As the product of litj^. 



(r)l 



'-'k(r) I 



appearing in (70) is only appreciably nonzero near the Fermi 
surface, the magnitude of the correlations is maximum there. 
This results in the ring-shape seen in both (a) and (b). 



For example if the detector in located in the x-y plane at 
a distance zq from the origin then the column integrated 
density-density correlation function is 



^0 

D{r±,r'^,t)^ / dzdz' D{r,r',t) 



(^\ 



6 _ 



— OO 

(70) 

where rj_ corresponds to spatial directions perpendicular 
to the column integration. As the p-wave breaks rota- 
tional symmetry, here chosen to be in the z direction, the 
observed column integrated density-density correlations 
depends on the relative orientation of the spontaneously 
broken symmetry direction and the detector. Figure 13 
shows the column integrated density-density correlation 
(weight function only) for a detector located in an x-y 
plane and an x-z plane for the p^ -\- ipy state with a gap 
Aq = 0.03ep>.. In Fig. 13b the p-wave nature can clearly 
be seen when the integrated column density is obtained 
by integrating along the direction that is perpendicular 
to the symmetry of the order parameter. 



^^'^^'^'^-Ltj 



6 ., 



) J(rmA + r'mA)l"k(r)l'kk(r)l'. 



(69) 

where S{r) is essentially a broadened delta function, the 
form of which depends on the specific details of the 
trapped system. Excluding the pre-factor, the weight of 
the "delta function" |u^.j.J ^|u,^,j. J ^ provides the informa- 
tion about the momentum correlations with k(r) = rm/t. 
Experimentally one does not measure the local density 
in 3-space, but instead the column integrated density. 



VII. CONCLUSIONS 

Cold atom experiments have demonstrated the capa- 
bility to study a remarkably simple many-body physics 
problem: That of two species of attractively interact- 
ing fermion as a function of the interatomic scattering 
length and relative densities of the two species. Despite 
the simplicity of this problem, the resulting phase dia- 
gram is quite rich, showing regions of phase separation, 
imbalanced superfluid, and normal Fermi liquid. 
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The question wc pursue here is, why aren't there more 
phases of imbalanced Fermi gases? Indeed, experiments 
on imbalanced Fermi gases observe the absence of any 
broken symmetry phases for a large range of parame- 
ters. For example, at unitarity, Ref. [38] finds that phase 
separation vanishes above P ~ 0.36. Arc there truly 
no brokcn-symmctry ground states of imbalanced Fermi 
gases over the range of polarization values 0.36 ^ P < I, 
or do other phases lurk at low temperatures in this 
strongly interacting system? 

This paper partially addresses such questions by 
proposing that, in the large imbalance region the true 
ground state is a p-wave superfluid of the spin-f fermions, 
with the order setting in below a temperature, Tc, plotted 
in Fig. 1. Important questions for futher work include 
obtaining more accurate estimates for Tc and studying 
the p-wave gap equation at low T, and determining the 
transition temperature for the onset of pairing for the 
spins-J, in the unitariy regime (and finding how this onset 
modifies the properties of the spin-t Cooper pairs) . An 
additional question concerns finding more experimental 
signatures of the onset of p-wave pairing, to help ascer- 



tain the validity of this scenario. 

From a general point of view, a natural question is 
whether other phases, such as higher-angular momentum 
superfluids or FFLO phases intervene in the T ^' limit 
of imbalanced Fermi gases. Since the p-wave interaction 
becomes repulsive for small imbalance (below the dashed 
line in Fig. 8), it is likely that other angular momentum 
channels can become dominant in this regime. Addition- 
ally, although the window of FFLO stability is extremely 
thin for the simplest FFLO-type state within mean-field 
theory, it is possible that generalized FFLO states, within 
theoretical approaches that go beyond mean-field theory, 
can lead to a wider regime of FFLO stability [21-23]. 
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